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Abstract 

Within the context of an effective potential formalism we calculate the contribution 
to Higgs boson masses in the next-to-minimal supersymmetric standard model from 
squark loops. We then supplement a previously performed renormalisation group 
analysis of the Higgs sector of this model with these results in order to determine the 
shift in the bound on the lightest CP-even Higgs boson mass as a result of squark 
effects. The improved bound on the lightest neutral CP-even Higgs boson mass, 
including squark contributions, is ruh < 146 (139, 149) GeV for = 90 (140, 190) 
GeV. For rut = 190 GeV squark effects contribute 23 GeV to the bound, with smaller 
contributions for small mt values. 



The next-to-minimal supersymmetric standard model [|l], |^ (NMSSM) provides an 



elegant extension of the minimal supersymmetric standard model (MSSM) which 
eliminates the /i-problem, the appearance of a mass-scale in the superpotential. The 
NMSSM differs from the MSSM by the presence of a gauge singlet field whose vacuum 
expectation value (vev) plays the role of the mass parameter /i in the MSSM, and is 
defined by the superpotential 



where the superfield = (tL,bL) contains the left-handed top and bottom quarks, 
and contains the charge conjugate of the right-handed top quark; Hi and H2 are the 
usual Higgs doublet superfields and is the Higgs gauge singlet superfield; the ellipsis 
represents terms whose relatively small couplings play no role in our analysis (in 
particular, we assume that the bottom quark Yukawa coupling, hb, satisfies hb <^ ht). 
The fields = {H^,Hi), H2 = {H21H2) and N develop vevs which may be 
assumed to be of the form 



where z/i, V2 and x are real, and y i/f + z/f = = 174 GeV. The low energy physical 
spectrum of the Higgs scalars consists of 3 CP-even neutral states, 2 CP-odd neutral 
states, and 2 charged scalars. A third CP-odd state is a Goldstone mode which 
becomes the longitudinal component of the while a further two charged scalars 
become those of the VT^'s. 

The tree-level bound on the lightest CP-even Higgs scalar mass in the NMSSM 




(1) 




< N >= X, 



(2) 
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ml<Ml + {X^u^ - M|) sin^ 2/3, 



(3) 



1 



where tan/3 = vijvx-, and radiative corrections to it 0, have been the subject of 
much discussion. In Ref. renormahsation group (RG) analyses were performed 
assuming either supersymmetric (SUSY) RG equations down to near the top quark 
mass, or the existence of only one light Higgs boson with a mass below the SUSY 
breaking scale. These analyses assumed degenerate squarks decoupling below their 
common mass scale Msusy- In Refs. within the effective potential formalism, 

squark effects have been considered in the NMSSM. However, no systematic attempt 
was made to study their effects on the bound in Eq. (j^), but rather one slice through 
parameter space was examined 0. Moreover, radiative corrections to the charged 
scalar mass-squared matrix were not presented, though these do not change the 
bound. 

Recently we performed a low energy renormahsation group analysis of the Higgs 
sector of the NMSSM and found that rrih < 146 GeV . This analysis considered the 
more general case of all the Higgs bosons having masses below the SUSY breaking 
scale (Msusy), and used non-SUSY RG equations below Msusv- However, in that 
analysis the squarks were assumed to be degenerate at Msusv = 1 TeV. Here we 
supplement this calculation with a computation of the stop and, where necessary, 
sbottom contributions to the CP-even, CP-odd and charged mass-squared matrices, 
from which we may determine modifications to the bound due to a general spectrum 
of squarks. We then attempt systematically to examine the squark corrected bound 
by surveying as much of the parameter space as practicable. 

Our previous analysis worked exclusively with Vo(/u), /i = 150 GeV being approx- 
imately the bound on the lightest CP-even state. From the superpotential in Eq. (||) 
we may obtain the pure Higgs part of Vq, where we have renormalised at the MS 
scale /i: 
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- X,\HlH,\' + X,\N\'\Hi? + X.Wm' 
+ Xj{N''^HiH2 + h.c.) + Xs\N\'' 

+ m\\Hi\^ + ml\H2\^ + ml\N\^ 

- mi{H^H2N + h.c.)-\m^{N^ + h.c.), (4) 
where Aj = Aj(/i) and mj = mi{fi). At Msusy the couphngs Aj and mj are given by: 

X, = X2 = {gl + gl)lA, X^ = {gl-gl)/A, 

A4 = A^ — (72/2, A5 = Ae = A^, A7 = — Afc, Ag = fc^, 
mi = mHi, 1712 = mH2, ms = mN, 1714 = XAx, = kAk. (5) 

(7i and g2 are the f/(l)y and SU{2)l gauge couphng constants, respectively, and m^/-, 
rriN, A\ and A^ are soft SUSY breaking parameters. We use the RG equations [0 to 
allow to us calculate the couplings renormalised at = 150 GeV given their values 
at Msusy- 

In order the calculate the shifts in the mass-squared matrices due to squarks, 
we employ the full one-loop effective potential. Much of our analysis of the squark 
spectrum is similar to that of the MSSM in Ref. whose notation we follow closely. 
The one-loop effective potential is given by 

ViW) = Vo(Q) + AVi(Q), 

where Vq{Q) is the tree-level potential at the arbitrary MS scale Str denotes 
the usual supertrace, and is the field-dependent mass-squared matrix. We shall 
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take Q = Msusy = 1 TeV. [] Vq is restricted to the pure Higgs part of the tree- 
level potential. Since our RG approach has already included the logarithmic effects 
of top quark and Higgs boson loops we exclude these particles from the supertrace. 
The logarithmic contributions to AVi{Q) from the top quark and Higgs bosons may 
be absorbed into Vo{Q) where Q = I TeV, to yield Vo(yu) where fi = 150 GeV, as 
discussed above. It only remains to calculate AVi{Q), where Q = 1 TeV, involving 
squark contributions. 

To calculate the shifts in the Higgs mass matrices as a result of squark effects 
we require the field-dependent mass-squared matrices of the squarks. Taking the 
contribution to the soft SUSY breaking potential involving squarks to be 

AVsoft = htAtiQH^i'^ + h.c.) + mjigp + ml\i^\^ + m\\¥\'^, (7) 

where V contains the charge conjugate of the right-handed bottom quark, and tildes 
denote the scalar components of the superfields, together with the superpotential in 
Eq. (P, leads to the field-dependent squark mass-squared matrix (ignoring contri- 
butions proportional to gauge couplings and hi,) 

( ml + hl\H^\'^ XhtNHl + htAtH^ -hlH^Ht \ 

2_ XhtNH'l + htAHl ml + hl{\Hl\'' + HtH:,) XhtN Ht - htAH^ 
-hfH^H^ XhtNH{ - htAtH^ + hlH^H^ 

^ m| y 

(8) 

in the basis {tL,t%hL,¥^, where a bar denotes complex conjugation. Manifestly, 
one eigenvalue is field-independent and may be discarded. Only the upper 2x2 
submatrix contributes to the CP-even and CP-odd mass-squared matrices, since the 
charged fields have zero vevs in order not to break QED, whereas the upper 3x3 
submatrix contributes to the charged mass-squared matrix. By differentiating AVi 

once with respect to the fields we may obtain the shifts in the minimisation conditions 

^This value of Q is of course totally arbitrary. The value of 1 TeV is chosen to be consistent with 
our previous analysis. The results are independent of Q to one loop order. 
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induced by squark loops, and twice will yield the shifts in the mass matrices. ^ 

In the basis {Hi, H2, N} we have the following mass-squared matrices, after ensur- 
ing that the full one-loop potential is correctly minimised. The couplings Aj are those 
obtained from the potential in Eq. (|^) renormalised at the scale fi = 150 GeV. For 
notational simplicity we drop the tildes. Moreover, we work in the basis of mass eigen- 
states {ti, t2, bi, 62} (though, since hb, the bottom quark Yukawa coupling, is taken to 
be zero, 62 never contributes to the mass-squared matrices), where the mass of the ti 
squark is m^j, and so on for the other squarks. The CP-even (scalar) mass-squared 
matrix is 

= + SM^, (9) 

where 



+ 



2Aiz/2 



2(A3 + A4)z/iI/2 



2(A3 + A4)z/iz/2 2A2Z/? 



2X5X^1 
tan i3[m4X — X-jX^ 



21^2 

2XqXI'2 



2A5XZ/1 

2X^x^2 
AXsx"^ — m^x 

[m^x — Xix^] —^[rriiX — 2X-jX^] 
— [rriiX — X^x'^] cot /3[m4X — X^x'^] —^[rriiX — 2A7X^] 
-^[m4a; - 2A7a;2] -^[rmx - 2Xjx^] ""^[rriix] 



(10) 



and 



/ a2 a2 a2 

' ^11 ^12 ^13 

A2 A2 A2 

^12 ^22 ^23 

^23 ^33 



V A2 
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The A'^j and are given by 



+ 



tan/3 -1 
-1 cot/3 



1^2 



Ui V1V2 



A2 



(11) 



A^ 

p 



A^ 
^11 



A^ 
^22 



A^ 
^33 



167r2 
3 



hl{Xx).At.f{ml,ml), 



r 4 2 m n2 Mt + Ax cot /3 
i^.-2.(Ax) . ^2^_^2^ 



9 {ml, ml), 



■/iV. fin ""'i""*^ ' 2A(A + Axcot/3) 



2'^t''2- 



^ -,4..2 



^SUSY 



+ 



ml - ml 



In^' 



m 



ti . 



2 f At{At + Xx cot PY 



mt 



mt 



9 (ml, ml). 



hy2-{Xiyi)\ 



' At + Xx cot /3' 
, ml - ml ^ 



9[mt^,mt^, 



^In fact, this is only approximately true due to Higgs self-energy corrections; these are expected 
to be small for the lightest Higgs bosons. 
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^12 



^13 



A^ 



87r2 
3 



At + Ax cot 



m| - ml 



In 



m 



t2 



+ 



At{At + Xx cot/5) 



mt 



m| - ml 



9 {ml, ml) , 



Stt 



At + Xx cot p\ ,0 2 \ 
9{mt,,rntJ 



ml - ml 



in 

^^2 



hl{Xx).{Xiyi).f{ml,ml), 



X 



12 5 



(12) 



where the functions / and g are defined by 

1 



fimt,,mt^] 
9 {ml, ml) 



ml - ml 



2 2 

m. In — ^-^^ — m. In — ^-^ — 



+ m 



t2 



? +m?Jln^ + 2('m? - 



mt 



(13) 



If the squarks are degenerate in mass, then equation Eq. (|) imphes that At + 
Xx cot f3 = 0. Furthermore, if m^^ = = Msusy, then Eq. (^) and Eq. (|T^ 
imply that the squark contribution to the CP-even mass-squared matrix vanishes. 
This is the hmit of our previous analysis. 

The CP-odd (pseudoscalar) mass-squared matrix is 



(14) 



where 



and 



tan P[m4X — X-jx"^] 
[m4X — Xjx^] 
^[m4X + 2XrX^] 



[m^x — Xix^] 
cot Plm^x — Xix"^] 
^[m^x + 2Xjx'^] 

tan (3 1 
1 cot (3 

V2 I'l 



f[miX + 2X7x'^] 
^[miX + 2Xjx^] 
3m5X + ^[m^x - ^X-jX^] 



V2 

X 

,.2 



A^ 

p 



Finally, the charged mass-squared matrix is 



tan (3 1 
1 cot/5 



{m4X — X-jX^ — A4Z/1Z/2 + A^ 



where 



A? 



IGtt^ 



E 



m„ In 



m„ 



d'^m^ 



^^■'-SUSY 



OH^dH; 



-)- I vevs 1 



(15) 



(16) 



(17) 
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and 



d^ml h^vliXxf. cot 13 hl{Xx).At 



2 ' 



dH^ dHt {ml - m|)(m?^ - mfj ml - mt^ 

d^ml _ h^ul{Xxf.cot(3 , hl{Xx).At 



+ 



dH^ dHt [ml - m^J(m| - m^^) ml - ml ' 

d^ml _ _ h^uliXxf.cotp 

dH{ dHt {ml - ml ) {ml - ml ) ' 



(19) 



The bound on the hghtest CP-even Higgs mass is a consequence of the fact that 
the minimum eigenvalue of M| is bounded by the minimum eigenvalue of the upper 
2x2 submatrix of M^. Using Eq. (|l^) to obtain the physical charged Higgs mass 
squared, m^, in terms of m4, 

m^x - Ayx^ = ]^{ml + Xiv"^) sin 2/5 - A^, (20) 

we may eliminate m^ from in favour of m^, and thus write the upper 2x2 
submatrix in the form 

Mf = M'2 + 5M'2, (21) 

where 

»^/2 / SAiz/i^ 2{X^ + Xi)uiU2 \ ^( tan/? \ ^( 2 , ^ 2-, ■ 

^ =[2{X, + X,)u,u, 2X,vl J + i -1 cot^ j2("^^ + ^^'^)^^"2^' 

(22) 



and 



It is clear from the form of the second term on the right-hand-side of Eq. (|23| ) that 
the factor — A^ does not change the bound at all, but merely serves, for a fixed 
tan/3, to shift the bound, when plotted as a function of mc, to the left or right. We 
may then drop this term, since its presence may be re-parametrised by a shift in the 
free parameter m^- In this way mf,^ is eliminated. 

From Eq. ( ^2]) and Eq. (^) it is a simple matter to determine the shift in the 
minimum eigenvalue of M'^ due to squark effects. Defining A = (M'^)ii, B = (M'^)i2 



and C = (M'^)22, the shift is given by 



{A-C){Ai,-Al,)+ABAl 
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(24) 



Because this shift is a one loop effect, the couphngs Aj in A, B and C may be evaluated 
at any renormalisation point, since the difference between a coupling evaluated at two 
different renormalisation points is also a one loop effect, thus giving an overall error 
at the two loop level, which we neglect in our approximation. 

In our previous analysis we used triviality limits on the couplings ht, A and k P] 
to determine the bound on the lightest CP-even Higgs boson mass not including 
general squark effects, that is, with 5M'^ = 0. Let the triviality limit on A for a given 
ht be Xmax- In table ^ we reproduce the values of ht and Xmax which generated the 
bound of our previous analysis — all other values, for a given top quark mass, 
resulted in smaller values of the bound. It transpires that A; = in all cases. In 
order to maximise the bound resulting from Eq. we take the limit rric = oo, 
thus eliminating parameter. As a function of rric the bound approaches its 

maximum asymptotically as rric oo. The approach is rapid, with rric ~ 200 GeV 
being a good approximation to = oo. Thus, the bound resulting from the purely 
formal procedure of taking the limit rric = oo is not unphysical [0. Had this not been 
the case, then certainly we would have an upper bound, but one perhaps not capable 
of realisation in an actual spectrum. 

We use the values of ht and Xmax in table |1] to calculate the shift in the bound on 
the lightest CP-even Higgs boson mass-squared given by Eq. (0). The shift in the 
bound on the mass is denoted by Srrih- This is a hideously complicated function of 
many parameters, including the squark masses rrit^ and rrit^, the soft SUSY breaking 

parameter At, the vev of the gauge singlet field x, and m^. As a function of rric the 

■^Taking k = gives rise to an axion in the physical spectrum. However, a small, non-zero value 
of k is sufBcicnt to give the would-be axion a mass of several tens of GeV. Such a value of k does 
not invalidate our calculations. 
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Table 1: Lightest CP-even Higgs mass bound in the NMSSM including general squark effects. In 
row 1 is the top mass, nit, in GeV; in rows 2 and 3 those values of ht^MsuSv) and \rnax{M su Sy) ^ 
which produce the bound in our previous analysis, this being in row 4, in GeV; in row 5 is the 
contribution to the bound resulting from a general squark mass spectrum rather than degenerate 
squarks with mass MsuSY, in GeV; in row 6 is the lightest CP-even Higgs mass bound including 
general squark effects, in GeV; in rows 7 and 8 are those value of rrit^ and At, in GeV, which 
generate the bound in row 6; nrit^ = 1 TeV. 
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940 
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1500 
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shift appears, numerically, to be maximised for rric small, typically between 100 GeV 
and 200 GeV, with the difference between the shifts at these two points being less 
than 1 GeV. Thus, we set = 200 GeV in the general squark contributions to the 
bound. Doing this enables us to retain the bound from our previous analysis derived 
from Eq. (^) in the limit rric = oo, and simply maximise the squark contributions 
separately and add them to our previous bound to yield the new bound rrih = -|- 
6mh. The x dependence is not too strong, with the ratio r = x/v typically taking 
values close to 10 in order to maximise the bound. Where this is not the case, the 
difference between the bound at its maximum, as a function of x, and that at r = 10 
is less than 1 GeV. Thus, we set r = 10, this giving a reliable indication of the 
maximum shift. We calculate the shift for various ranges of squark masses and a 
range of values of At, and record the maximum value obtained. The squark masses 
are allowed to vary between 20 GeV and 1000 GeV in steps of 10 GeV (1000 GeV is 
the upper limit since we take Msusv = 1000 GeV.) Given this restriction, the value 
of At which maximises the general squark contribution to the bound never exceeds 2 
TeV. We impose the constraint 2hti'2{At + Ax cot (3) < |m^^ — mf^ \ which follows from 
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the form of the squark mass matrix in Eq. (|^). 

In table |l| and figure 1 we show the bound on the hghtest CP-even Higgs mass 
ruh including general squark effects (and m° for comparison). Table also shows the 
values of mt^ and At used to generate the bound. As rrit increases, monotonically 
decreases, while At monotonically increases. The other squark mass, rrit^, takes the 
value 1 TeV over the whole range of top quark masses — that one squark mass takes 
on its maximum permitted value to maximise the general squark contribution to the 
bound can be seen analytically. It will be noticed that these squark masses and At 
are somewhat larger than typically expected from Grand Unification scenarios. How- 
ever, we feel it wise not to restrict ourselves too excessively to particular prejudices 
regarding physics as yet unknown (though, of course, the notion of triviality does 
require the assumption of a SUSY desert up to the unification scale). 

Our calculations suggest that the universal upper bound on the lightest CP-even 
Higgs mass is lifted from 145 GeV to 149 GeV, a mere 4 GeV difference. Table |l] 
reveals that for larger values of rrit the shifts may be up to 23 GeV, but for small 
rut the squark effects are small, as expected. We emphasise that we have not re- 
maximised the complete one loop corrected bound, but rather made the assumption 
that the values of ht and Xmax which generated the bound of our previous analysis are 
not significantly modified by the inclusion of a general spectrum of squark masses. 

To conclude, we have supplemented our previous RG analysis of the Higgs sector 
of the NMSSM with a calculation of general squark effects, and determined how they 
modify the bound on the mass of the lightest CP-even state. Our calculations indicate 
that the effects of a general squark spectrum can be very significant for large top quark 
masses. One reason why, for large top quark masses, squark effects are important is 
that there exist finite, one-loop diagrams with vertices containing the uncontrolled 
soft SUSY breaking parameter At. These diagrams can give large contributions to 
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the Higgs boson masses. We close by indicating that in the NMSSM there are similar 
diagrams involving loops of Higgs bosons, again with vertices containing soft SUSY 
breaking parameters, this time Ax and A^. The existence of these diagrams is directly 
due to the gauge singlet field A^; they do not occur in the MSSM. We see no reason in 
principle why these new diagrams should not give rise to similarly large effects, except 
for obvious factors of 3 due to colour. We are currently in the process of estimating 
these effects [|T0| . 
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Figure Caption 

Figure 1: The bound on the Hghtest CP-even Higgs boson mass in the NMSSM. The 
upper line is the bound ruh including a general spectrum of squark masses. The lower 
line is the bound of our previous analysis, from assuming degenerate squarks with 

a common mass Msusy- 
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